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Abstract 

The dynamic impact response of giant buckyball C 72u is investigated by using molecular dynamics simulations. The 
non-recoverable deformation of C 72 o makes it an ideal candidate for high-performance energy absorption. Firstly, 
mechanical behaviors under dynamic impact and low-speed crushing are simulated and modeled, which clarifies 
the buckling-related energy absorption mechanism. One-dimensional C 72u arrays (both vertical and horizontal 
alignments) are studied at various impact speeds, which show that the energy absorption ability is dominated by 
the impact energy per buckyball and less sensitive to the number and arrangement direction of buckyballs. 
Three-dimensional stacking of buckyballs in simple cubic, body-centered cubic, hexagonal, and face-centered cubic 
forms are investigated. Stacking form with higher occupation density yields higher energy absorption. The present 
study may shed lights on employing C 720 assembly as an advanced energy absorption system against low-speed 
impacts. 
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Background 

Absorption of external impact energy has long been a 
research topic with the pressing need from civil [1,2] to 
military needs [3,4]. In particular, effective absorption of 
mechanical energy at low-impact speed, i.e., below 100 
m/s is of great interest [5,6]. As one of the major 
branches of fullerene family, the carbon nanotube 
(CNT) has demonstrated an outstanding mechanical en- 
ergy dissipation ability through water- filled CNT [7], 
CNT forest and bundle [7], CNT/epoxy nanocomposites 
[8], CNT immersed in nonaqueous liquid [9], intercalat- 
ing vertical alignment with aligned existing layered com- 
pounds [10], and sponge-like material containing self- 
assembled interconnected CNT skeletons [11], among 
others. The advantage lies within the CNTs' intriguing 
mechanical properties, i.e., ultra-strong (Youngs modu- 
lus of 0.9 to 5.5 TPa [12-14] and tensile strength of 60 
GPa [12]) and ultra-light, as well as the tube structure 
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which buckles upon external loadings [15]. Both theoret- 
ical modeling [16-18] and experiments [19-21] have pro- 
posed that the energy dissipation density of CNTs could 
be on the order of 200 J/cm 3 , about 1-2 order of magni- 
tudes over traditional engineering material [1]. 

Naturally, another branch of fullerene family with a 
spherical shape, i.e., the buckyball, also possesses excel- 
lent mechanical properties similar to CNTs. Man et al. 
[22] examined a C 60 in collision with a graphite surface 
and found that the C 60 would first deform into a disk- 
like structure and then recover to its original shape. It is 
also known that C 60 has a decent damping ability by 
transferring impact energy to internal energy [23,24]. 
This large deformation ability under compressive strain 
of C 60 was also verified by Kaur et al. [25]. For higher 
impact energy, Zhang [26] employed C 60 /C 32 o to collide 
with mono/double layer graphene, and the penetration 
of graphene and the dissociation of buckyball were 
observed. Furthermore, Wang and Lee [27] observed a 
novel phenomenon of heat wave propagation driven by 
impact loading between C 60 and graphene which was re- 
sponsible for the mechanical deformation of the bucky- 
ball. Meanwhile, giant buckyballs, such as C 72 o> have 
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smaller system rigidity as well as non-recoverable 
morphology upon impact, and thus they are expected to 
have higher capabilities for energy dissipation [28] . How- 
ever, to the best knowledge of the authors, currently, 
only few studies about the mechanical behavior of giant 
buckyball are available [29-31]. 

To understand the mechanical behavior of C 72 o and 
investigate its energy absorption potential in this paper, 
the dynamic response of C 72 o is studied at various im- 
pact speeds below 100 m/s by employing molecular dy- 
namics (MD) simulations. Firstly, the buckling behaviors 
under both low-speed crushing and impact are dis- 
cussed and described using classical thin shell models. 
Next, 1-D alignment of C 72 o system is investigated to 
identify the influence of packing of the buckyball on 
unit energy absorption. Finally, 3-D stacking of C 72 o sys- 
tem is considered, where four types of packing forms 
are introduced and the relationship between unit energy 
absorption and stacking density are elucidated by an 
empirical model. 

Methods 

Computational model and method 

The C 720 is a spherical buckyball with diameter of 2.708 
nm (where the van der Waals equilibrium distance is con- 
sidered), volume of 7.35 nm 3 , and mass of 1.45 x 10" 20 g. 
C 720 with varying numbers and packing directions (both 
vertical and horizontal) are selected in this study. Com- 
putational cells from single buckyball to 3-D buckyball 
stacking system are illustrated in selected examples in 
Figure 1. In the scenario of the impact, the buckyball 
system subjects to the impact of a top rigid plate with in- 
cident energy £i mpa ctop an d the initial impact speed is 
below 100 m/s; in the scenario of crushing, the top rigid 



plate compresses the buckyball system at a constant 
speed below 100 m/s. The bottom plate, which is rigid 
and fixed, serves as a receiver, and the force history it 
experiences could indicate the energy mitigation cap- 
ability of the protective buckyball system. The buckyball 
is not allowed to slip with respect to the impactor and 
receiver plates. Both the impactor and receiver plates 
are composed of carbon atoms. The masses of the atoms 
are varied in the following simulation to set various 
loading conditions (varying impactor mass), while the 
interactions between the plates and buckyballs remain 
as carbon-carbon interaction. 

MD simulation is performed based on large-scale 
atomic/molecular massively parallel simulator platform 
with the micro-canonical ensembles (NVE) [32] after 
equilibration. A pairwise Lennard-Jones (L-J) potential 
term is added to the buckyball potential to account for 
the steric and van der Waals carbon-carbon interaction 



Ufa) =4e cc 



Q"cc 



12 



free 
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where e C c is the depth of the potential well between 
carbon-carbon atoms, o C c is the finite distance where 
the carbon-carbon potential is zero, r t j is the distance be- 
tween the two carbon atoms. Here, L-J parameters for 
the carbon atoms of the buckyball Oqq = 3.47 'A and 
e C c = 0.27647 kj/mol as used in the original paramet- 
rization of Girifalco [33] and van der Waals interaction 
govern in the plate-buckyball interaction. A time inte- 
gration step of 1 fs is used, and periodical boundary con- 
ditions are applied in the x,y plane to eliminated the 
boundary effect. 




Receiver 



(a) 0-D Single C720 buckyball energy absorption system; 

(b) 1-D chain-like C720 buckyball energy absorption system with 3 buckyballs 
included 

(c) 3-D C720 buckyball energy absorption system in Square Cubic (SC) stacking 
form with 2 by 2 by 3 buckyballs included (x-z plane side-view) 

Figure 1 Various alignments of buckyball system as a protector. 
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Single buckyball mechanical behavior 
Atomistic simulation result 

The distinctive mechanical behavior of a single buckyball 
should underpin the overall energy absorption ability of a 
buckyball assembly. The force F and displacement W are 
normalized as FR/Eh 3 and W/D, respectively, where R, h, 
D, and E are the radius, effective thickness, diameter, and 
effective Youngs modulus of the buckyball, respectively. 
Considering that bending is involved during the buckyball 
compression, h = 0.66 nm and E = 5 TPa [34,35]. Here a 
crushing speed at 0.01 m/s is employed to mimic quasi- 
static loading, because the normalized force-displacement 
curves are verified to be the same at various loading rates 
from 0.1 to 0.001 m/s in trial simulations. The force- 
displacement response under both quasi-static and a rep- 
resentative dynamic impact loading (with impact speed of 
50 m/s and energy of 1.83 eV) are studied, as shown in 
Figure 2. Two obvious force-drops could be observed in 
low-speed crushing, while only one prominent force-drop 
exists in dynamic loading which is related to the less- 
evident snap-through deformation shape. 

The entire compression process could be divided into four 
phases according to the FR/Eh 3 ~ W/D curve, i.e., buckling 
(W/D < 10%), post-buckling (10% < W/D< 30%), densifica- 
tion (30% < W/D < 40%), and inverted-cap-forming phase 
(W/D > 40%). Upon the ricochet of the plate, the deform- 
ation remains as a bowl shape with great volume shrinkage. 
The stabilization of such a buckled morphology is owing to 
a lower system potential energy in the buckled configuration 
due to van der Waals interaction; similar energy dissipation 
mechanism in CNT network is also revealed by [36]. 

The derivative of curve undergoes a sudden change 
at the same W/D value but in two completely different 




Phase I Phase II Phase III Phase IV 
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Normalized displacement, W/D 

Figure 2 Normalized force displacement curves at both low- 
speed crushing and impact loading. The entire process from the 
beginning of loading to the bowl-forming morphology can be 
divided into four phases. Morphologies of C 720 are shown at the 
corresponding normalized displacements. 



loading rates, suggesting that the sudden force-drop 
points are highly dependent on the buckyball deformation 
rather than the loading rate. And theoretical insights may 
be obtained from the four-phase deformation. 

Phenomenological mechanical models 

Note that due to the property of FR/Eh 3 ~ W/D curve, 
among the phases of compression process, those with 
significant reduction of force (Figure 2) are relatively un- 
important for energy absorption and not included in the 
modeling effort. A three-phase model for low-speed 
crushing and a two-phase model for impact loading are 
proposed separately in the following sections. 

Three-phase model for low-speed crushing (quasi-static 
loading) 

(1) Phase I. Buckling phase 

In the range of small deformation in the beginning of 
compression, the model describing thin-shell deformation 
under a point force is applicable [37,38]. Considering a 
buckyball with wall thickness h = 0.066 nm compressed 
by F with deformation of W (with the subscript number 
denoting the phase number sketched in Figure 3), the 
force-deflection relation should be expressed as [39] 

F 1 = — W 1 (0<Wi<W r M ) (2) 
Rc 

where the bending stiffness G = Ehc 2 ; the reduced wall 
thickness c = hj \/l2(l — v 2 ) and v is the Poissons ratio. 
The linear deformation behavior continues until it reaches 
the critical normalized strain W^i* Experimental results 
for bulk thin spherical shell show that the transition from 
the flattened to the buckled configuration occurs at a de- 
formation close to twice the thickness of the shell [40]; 
while Wbi here is about 4/z, indicating a larger buckling 
strain in nanoscale structure. 

The nanostructure has higher resistance to buckle 
than its continuum counterpart and based on the fitting 
of MD simulation, a coefficient / ~ 2.95 should be 
expanded to Equation 2 as 

RC 

Fi=^W r 1 -(f*)(0< Wi< W bl ) (3) 
Rc 

It is reminded that this equation is only valid for C 72 o 
under low-speed (or quasi-static) crushing. 

(2) Phase II. Post-buckling phase 

As the compression continues, buckyball continues to 
deform. Once the compressive strain reaches W\>i, the 
flattened area becomes unstable and within a small 
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Low-speed crushing deformation mode 
Phase II 
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Figure 3 Illustration of deformation phases during compression for C 72 o. Dynamic loading and low-speed crushing share the same 
morphologies in phase I while they are different in phase II. Analytical models are based on the phases indicated above and below the dash line 
for low-speed crushing and impact loading, respectively. 



region, the buckyball snaps through to a new configur- 
ation in order to minimize the strain energy of the de- 
formation, shown in Figure 3. The ratio between the 
diameter and thickness of buckyball is quite large, i.e., 
D/h ~ 36.5, and only a small portion of volume is 
involved thus the stretching energy is of secondary order 
contribution to the total strain energy. Hubbard and 
Stronge [41] developed a model to describe the post- 
buckling behavior of a thin spherical shell under com- 
pression based on Steele's [42] model 



F 2 = ^^ 1 -^(W bl <W 2 lW b2 ) 



(4) 



(3) Phase III. Densification phase 

When the compression goes further, the crushing dis- 
placement eventually becomes much larger than the 
thickness and thus the force-displacement relation 
becomes nonlinear [42]. The buckled buckyball is densi- 
fied during this process. A phenomenological nonlinear 
spring-like behavior could be fitted as 



F 3 = yWg, 



(6) 



where y is a coefficient and n is fitted as n ~ 1.16. Con- 
sidering the relationship [41,42] 



where K = (|) (^) \/3(l — v 2 ). This nonlinear deform- 
ation behavior extends until it reaches the densification 
critical normalized strain W h2 . The value of W h2 could be 
fitted from the simulation data for C 720 where W\> 2 ~ 

The first force-drop phenomenon is obvious once the 
buckling occurs where the loading drops to nearly zero. 
Therefore, by applying the boundary condition of F 2 (W 2 ) ~ 
0, Equation 4 maybe further modified as 



h I 2 



and 
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(8) 



we may come to the equation 



l2hl6G / 
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Thus, by considering the continuity of two curves in 
adjacent phases, we may rewrite Equation 9 as 



fl6Gh ( 2 



-F 2 (W b2 )y(f )(W b2 < W 3 ) 
(10) 



Therefore, Equations 3, 5, and 10 together serve as the 
normalized force-displacement model which may be used 
to describe the mechanical behavior of the buckyball under 
quasi-static loading condition from small to large deformation. 

Figure 4 shows the simulation data at low-speed crush- 
ing compared with the model calculation. A good agree- 
ment between two results is observed which validates the 
effectiveness of the model. 

Two-phase model for impact 

The mechanical behaviors of buckyball during the first 
phase at both low-speed crushing and impact loadings 
are similar. Thus, Equation 2 is still valid in phase I with 
a different / ~ 4.30. The characteristic buckling time, 
the time it takes from contact to buckle, is on the order 
of t « 10" 1 ~ 10° ns ~ T~ 2.5R/C! » 5.71 x 10" 5 ns, where 
p is the density of C 72 o and C\ = y/E/p . It is much 
longer than the wave traveling time; thus, the enhance- 
ment of / should be caused by the inertia effect [43] . 

As indicated before, the buckyball behaves differently 
during the post-buckling phase if it is loaded dynamic- 
ally, i.e., no obvious snap through would be observed 
at the buckling point such that the thin spherical 
structure is able to sustain load by bending its wall. 
Therefore, a simple shell bending model is employed 
here to describe its behavior as shown in Figure 3; the 
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Figure 4 Comparison between computational results and 
analytical model at low-speed crushing of 0.01 m/s. 
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Figure 5 Comparison between computational results and 
analytical model. Comparison between computational results and 
analytical model at various impact speeds from 40 to 90 m/s. 



top and bottom flattened wall with length of L experi- 
ences little stretching strain, whereas the side wall 
bends with finite deformation, governing the total sys- 
tem strain energy 



-^system 
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(ii) 



and M is the 



where the bending rigidity EI 

bending moment. A denotes the integration area. The 
h is the 'enlarged' thickness, the result of smaller 
snap-through phenomenon. Here, h ~ 1.40/z via data 
fitting. Substituting geometrical constraints and taking 
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Figure 6 Characteristic normalized force-displacement curve of 
1-D system with vertically lined C 72 o buckyballs. The 

characteristic normalized force-displacement curve of 1-D system 
with five vertically lined C 720 s at impact speed of 50 m/s. 
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Normalized displacement, W/D 

Figure 7 Characteristic normalized force-displacement curve of 
1-D buckyball system with various numbers of horizontally 
lined C 720 buckyballs. The characteristic normalized force- 
displacement curve of 1-D buckyball system with various numbers 
of horizontally lined C 720 s at impact speed of 50 m/s. 
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Figure 9 UME and maximum contact force at constant impact 
speed (50 m/s) with various impact masses. UME and maximum 
contact force at constant impact speed (50 m/s) with various impact 
mass (from 8.7 x 10~ 19 to 7.1 x 10~ 17 g), and constant impact mass 
(2.8 x 10" 18 g) with various impact speeds (from 10 to 90 m/s), for 
five-buckyball systems. 



the derivative, the force-displacement relation becomes 
(for C 72 o under 100 m/s impact) 



-EI ; 

4 (R- W 2 /2) 



jt 2 -jtR (W bl < W 2 <W b2 ) 



(12) 



Therefore, Equations 3 and 12 together provide a 
model to describe the mechanical behavior of the bucky- 
ball under dynamic loadings. 

When the impact speed is varied, the corresponding 
force is modified by a factor a owing to strain rate effect 
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Figure 8 UME and UVE values of both vertical and horizontal 
buckyball systems with various buckyball numbers. UME and 

UVE values of both vertical and horizontal buckyball systems with 
various buckyball numbers at impact speed of 50 m/s. 



[44-46]. With the subscript representing the impact speed 
(in units of m/s), the correction factor c = [a 40 , a 50 , a 60 , a 70 , 
ttso> ago] = [0-83, 1.00, 1.12, 1.14, 1.17]. Figure 5 illustrates 
the comparison between atomistic simulation and model 
(for impact speeds of 40 to 90 m/s), with good agreements. 

Results and discussion 

Buckyball assembly 

In practice, buckyballs need to be assembled (shown in 
Figure 1) so as to protect materials/devices. Various stack- 
ing arrays are investigated as follows. 
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Figure 10 Normalized force-displacement curves for SC BCC 
FCC and HCP packing of C 720 . Typical normalized force- 
displacement curves for SC, BCC, FCC and HCP packing of C 72 o at 
impact speed of 50 m/s, and the impact energy per buckyball is 
1 .83 eV. 
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1-D alignment buckyball system 

The C 72 o can be arranged both vertically and horizon- 
tally in a 1-D chain-like alignment. Figure 6 shows the 
mechanical behavior of a five-buckyball array subjecting 
to a rigid plate impact with impact energy and speed of 
9.16 eV and 50 m/s respectively. Progressive buckling 
and bowl-shape forming behavior takes the full advan- 
tage of single buckyball energy absorption ability one by 
one and controls the force on the receiver within a rela- 
tively low value during first section of deformation 
(within W/D < 1.5) which provides cushion protections. 

Another 1-D arrangement direction is normal to a plate 
impact. Unlike the progressive buckling behavior in the 
vertical system, all buckyballs buckle simultaneously in 



the horizontal array. Figure 7 shows the scenario with im- 
pact energy of 1.83 eV per buckyball and impact speed of 
50 m/s, where the total reaction force scales with the 
number of buckyballs. Systems with different buckyball 
numbers show almost uniform deformation characteristics 
of individual buckyballs. 

The energy absorption per unit mass (UME, J/g) and 
unit volume (UVE, J/cm 3 ) are given in Figure 8, which 
shows that the UME and UVE are almost invariant re- 
gardless of buckyball number or arrangement. In Figure 8 
the impact energy per buckyball is fixed as 1.83 eV; if the 
impact energy or speed changes, the value of UME or 
UVE alters; however, the result is still insensitive to bucky- 
ball number or arrangement. The major responsible 
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Figure 11 UME and UVE values of SC, BCC, FCC and HCP packing of C 720 at impact speeds. UME and UVE values of SC, BCC, FCC, and 
HCP packing of C 72 o at impact speeds from 10 m/s to 90 m/s. Fitting surfaces based on the empirical equations are also compared with the 
simulation, (a) UME values of various packing forms of C 720 at impact various impact speeds, (b) UVE values of various packing forms of C 720 at 
impact various impact speeds. 
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reason is that the energy absorption ability of the system 
stems from the non-recoverable deformation of individual 
buckyball which is almost uniform. 

By fixing either the impact speed or mass and varying 
the other parameter, the impact energy per buckyball 
can be varied. It imposes a nonlinear influence on the 
UME and the maximum force on the receiver, as shown 
in Figure 9 for the vertical alignment of five-buckyball 
system. No matter how the impact speed or mass varies, 
it is the impact energy per buckyball that dominates the 
values of UME and maximum transmitted force. 



3-D stacking buckyball system 

The packing density of a 3-D stacking system can be dif- 
ferent than that of the 1-D system, and thus the perform- 
ance is expected to vary. Four types of 3-D stacking forms 
are investigated, i.e., simple cubic (SC), body-centered 
cubic (BCC), face-centered cubic (FCC) (a basic crystal 
structure of buckyball [47]), and hexagonal-closed packing 
(HCP). The occupation density // 5C = 7r/6 « 0.52, rj BCC = 
71^/3/8-0.68, rj ¥CC = rj HCP = tt/3v^*0.74 [48] for SC, 
BCC, FCC, and HCP, respectively. Convergence study 
indicates that the profiles of force-displacement curves as 
well as the energy absorption rate at increasing buckyball 
numbers at one computational cell keep the same. In this 
case, a fundamental unit, such as containing 2x2x3 
buckyballs for SC arrangement is shown in Figure lc. 

Figure 10 illustrates the normalized force-displacement 
curves for SC, BCC, FCC, and HCP units under the same 
impact energy per buckyball (1.83 eV). As expected, the 
mechanical behaviors of FCC and HCP are similar, while 
the BCC and SC units (with lower rj) have more space for 
system to comply and hence the impact force is smaller 
yet the displacement is larger. Consequently, FCC and 
HCP have the same energy absorption ability and that of 
BCC and SC are inferior. 

Energy absorption performances of the three basic units 
are studied at various impact speeds, i.e., from 10 to 90 
m/s while the impact mass is kept a constant, as shown in 
Figure 11. With the impact speed increases, more mech- 
anical energy is absorbed; but the increasing trend 
becomes slighter at higher impact speed when the bucky- 
ball system reaches its mitigation limit. The improvement 
is greater in terms of UVE than UME with higher rj. 

By normalizing the UME and UVE as 77 m = UME/ 
(£impactor/^) and 77 v = UVE/^impactor/Koiume) where 
Koiume is the volume of the buckyball and impact speed 
as V = v/y r B/p where B = 34 GPa [49] is the bulk 
modulus of graphite. An empirical equation could be fit- 
ted as 

n m =Arj(BV c +DV) (13) 



where A = 5.50, B = -0.25, C = 0.21, and D = 25.0 with 
fitting correlation coefficient of 0.96 and 

n v =Arj(BV c +DV), (14) 

where A = 0.46, B = -1.94, C = 0.21, and D = 187.9 with 
fitting correlation coefficient of 0.96. These equations 
are valid for low-speed impact speed (below 100 m/s) on 
stacked C 72 o buckyballs. When the impact speed is fixed, 
the unit energy absorption linearly increases with the oc- 
cupation density; under a particular spatial arrangement, 
the energy absorption ability increases nonlinearly with 
the impact speed. 

Conclusions 

C720 as a representative giant buckyball has the distinct- 
ive property of non-recovery deformation after crushing 
or impact, which makes it capable of absorbing a large 
amount of energy. The mechanical behaviors of a single 
C720 under quasi-static (low-speed crushing) and dy- 
namic impact are investigated via MD simulation and 
analytical modeling. By understanding the mechanism of 
mechanical behavior of individual C 72 o, the energy ab- 
sorption ability of a 1-D array of buckyball system is 
studied. It is found that regardless of the direction of 
alignment and number of buckyballs, the unit energy ab- 
sorption density is almost the same for low-speed im- 
pact. In addition, different 3-D stacking at various 
impact speeds and stacking forms are investigated. Ex- 
plicit empirical models are suggested where packing 
density and impact speed may pose a positive effect on 
the unit energy absorption. This study may shed lights 
on the buckyball dynamic mechanical behavior and its 
application in energy absorption devices and inspire the 
related experimental work. 
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